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Abstract 

A good fit of the fermion masses and mixings has been found in the 
minimal renormalizable supersymmetric SO(10). This solution needs a 
strongly split supersymmetry breaking scenario with gauginos and higgsi- 
nos around 100 TeV, sfermions close to 10 GeV and a low GUT scale 
of around 6 x 10 15 GeV. We predict fast proton decays through SO(10) 
type of d = 6 operators and the leptonic mixing angle sin 6*13 0.1. 

1 Introduction 

Assuming renormalizability and no other symmetries, the supersymmetric 
SO(10) with three generations of matter 16/? and the Higgs representations 
of 210//, 126//, 126// and 10// [H[2] turns out to have only 26 free parameters 
(on top of the usual soft supersymmetry breaking terms) [3] and thus can be 
considered the minimal prototype grand unified theory (GUT) model. 

Several efforts have been employed in trying to fit the fermion masses in 
this minimal SO(10) [2 El El El EJ QUI EH E2] . It has been finally found out 
that due to constraints from the Yukawa and Higgs sectors either the neutrino 
mass scale came out too small or the gauge coupling constants entered the 
non-perturbative regime [131 Ell EI3 EE]- 

Two main objections can be raised to most of the above works. 

First, the gauge and Yukawa couplings at the GUT scale were obtained 
from the assumption of a desert from the low energy supersymmetry (SUSY) 
scale to the GUT scale. This is not always true, since, depending on the 
various Higgs couplings of the model, there could be states lying one or two 
orders of magnitude below the GUT scale. In a scenario with small SO(10) 
representations [IHEE1I19, 20J this would always be just a small perturbative 
correction, but in the minimal SO(10) the large number of light remnants of 
the 210//, 126// or 126// representations could make the difference. Although 
in some cases it has been found that the relevant threshold corrections are 
negligible [21], this has not been studied systematically. 



Second, the unknown values of the soft terms in the minimal supersym- 
metric standard model (MSSM) make any prediction of the fermion masses 
and mixings impossible. Due to (in principle) large finite threshold correc- 
tions at the low supersymmetry breaking scale [221 [22] (the possible impact 
of these effects for renormalizable SO(10) has been recently studied in [24J) 
it is in general difficult to find any fitting reliable. Even if a good fit were 
found, unknown soft terms would change it. So to have a really predictive 
model one would need to know how supersymmetry is broken and mediated. 
Of course there are many such models on the market, but their construction 
is usually completely orthogonal to the GUT one is considering (an opposite 
attempt to construct such a supersymmetry breaking model in a GUT has 
been recently given in [25J). Also, due to the same reason, d — 5 proton 
decay is not predicted by the theory, although it will typically be danger- 
ously large. In short, unless the finite threshold corrections are known to be 
negligible for some reason, there is little hope to have a grand unified theory 
of fermion masses and proton decay. 

Here we will show that a consistent treatment of the first issue auto- 
matically pushes to the split supersymmetry (26J, E7J EH] scenario, due to 
which the second problem is automatically solved: threshold corrections to 
the fermion masses are suppressed by powers of the large sfermion masses. 
Obviously, as usually in the split SUSY scenario, there are no d = 5 proton 
decay operators, as well as no dangerous flavour changing neutral currents in 
the SM. The only solution is possible however for relatively small values of 
the GUT scale, very close to the experimentally allowed value. If this model 
turns out to be correct, the next generation of proton decay searches must 
be unavoidably successful Q Not only this: the presence of d = 6 operators 
only and the known fermion flavour mixing matrices obtained from the fit 
make proton decay well known in all its channels. The model thus connects 
proton decay with neutrino masses, one of the ultimate goals of all grand 
unified theories. 

The paper is written as follows: in section [2] we will shortly review the 
minimal renormalizable supersymmetric SO(10) model, as well as set up the 
notation to be used later. In section [3] we will first summarize the (well 
known) main problem in the fitting so far: the small neutrino overall scale. 
We will then propose a general strategy of how to remedy that problem in 

1 This conclusion may change at the 2-loop order: after all, the fit determines only 
log Mgut, while the proton decay lifetime is proportional to Mq ut . 
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a phenomenologically acceptable way. First, from gauge coupling unification 
and d — 6 proton decay constraints the allowed parameter space will be de- 
duced. Second, we will scan that parameter space by performing a general 
fit of all fermion parameters except the neutrino overall scale. In this way we 
will see how close this scale can be to the one in agreement with experiments 
(in all interesting fits the neutrino spectrum will turn out to be normally 
hierarchical) for acceptably small x 2 of the general fit. We will conclude that 
a factor of about 5 will be missing in order to get the correct neutrino mass 
scale in a scenario with low energy supersymmetry, thus confirming the con- 
clusions of [131 IH US! US] • This will push us to consider in section [i] a scenario 
with a different sparticle spectrum: instead of low energy supersymmetry we 
will allow the sfermion masses to differ from the gaugino and higgsino ones. 
This time our strategy will yield a successful numerical fit, predicting ob- 
servable proton decay rates, the supersymmetry breaking scale to be around 
10 13 -10 14 GeV with gauginos and higgsinos in the 100 TeV region and lep- 
tonic mixing angle sin# 13 « 0.1. A theoretically interesting outcome of this 
analysis is that neutrino masses get the main contribution from the type II 
seesaw mechanism. Finally we will outline the main results in section [5] and 
describe some further work to be done. 



2 The minimal SO(10) 

The minimal supersymmetric SO(10) contains on top of the usual three gen- 
erations of spinorial 16 p also the 210 h, 126 h, 126h and 10 h Higgs represen- 
tations. The renormalizable superpotential W is 

W = 16 F (Y 10 10„ + Y m T26„)16 F 

+ -210| + -2101, + —126^126^ + ^126^210^126^ 

+m^l0 H 2 + ^210^10^(al26 H + aT26H), (1) 

where Y\q and Y^ are the two complex symmetric Yukawa matrices of the 
theory. 

The above Lagrangian will be used for the determination of the mass 
spectrum needed for the renormalization group equation running of gauge 
and Yukawa coupling constants on one side and for the derivation of the 
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relevant mass matrices needed for the fitting of the light fermion masses and 
mixing parameters at the GUT scale on the other side. 

For the first one of these two topics, i.e., the mass spectrum, we refer the 
reader to the existent literature [29], EH EE] (for very recent reanalysis see 
[3H E2] ) - The upshot of those studies is that after a required fine-tuning the 
Higgs sector has only eight real parameters: 

m, a, a, |A|, \rj\, 4> = arg(A) = — arg(^), x = Re(x) + ilm(x). (2) 



Here we closely follow the notation advocated in [H]. Once these parameters 
are given the mass spectrum of the Higgs fields is completely determined and 
can be found numerically. (We set = for simplicity in what follows.) 
These parameters also specify all the vacuum expectation values (VEVs) of 
the theory and hence, besides Yukawa couplings Y w and Yj^q, directly affect 
the predictions for the fermion masses and their mixing parameters. 

Let us now set up the notation regarding the fermion mass matrices fol- 
lowing [14] . All we will need are two matrix sum-rules: 

N 

M u = -±taaP[(l + £(z)W d -Z(x)M e ], (3) 

v sin 2 /? /lAliV 2 

M n = % a \ rTTT m ^ J x )+™nfii(x , 4 

m cos p V \r)\ N d 

where N u and N d as well as all the ratios of polinomials of x, i.e., //, fn 
and £, are specified in Ref. [H]. tan/3, another parameter of the model, 
is the ratio of up-type and down-type VEVs of the MSSM-like Higgs fields 
and v(— 174 GeV at the Mz scale) is the scale of SU(2) breaking, mj = 
M e (M d - M e )~ x M e - 6£M e + % 2 {M d - M e ) (m n = M d - M e ) is proportional 
to the type I (type II) seesaw contribution to the light neutrino mass matrix. 
We remind the reader that type I seesaw is mediated by the right-handed 
neutrinos [33l [Ml ESI [36l ETj, while the type II seesaw comes from the triplet 
VEV [3H1 HSJ HOl HI] . M u (d) and M e („) are 3x3 complex symmetric mass 
matrices of up (down) quarks and charged (neutral) leptons, respectively. 
Eqs. @ and Q are crucial to establish whether a viable description of known 
fermion masses, mixing angles and phases at the GUT scale is possible. We 
will accordingly use them whenever we perform numerical fits. 
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3 The neutrino mass scale 



The crucial question is whether this model can generate heavy enough light 
neutrino masses without being in conflict with other phenomenological con- 
straints. 

3.1 Maximization of the neutrino mass scale 

This issue is especially easy to illustrate in the context of a dominant type 
II seesaw mechanism. We accordingly focus on the magnitude of the factor 
that multiplies the mismatch between the down-quark and charged lepton 
mass matrices, i.e., Ma — M e , and thus controls the strength of the type II 
contribution 



We want it to be as large as possible and still phenomenologically viable. We 
will comment on the type I contribution in concrete examples to show that 
due to gauge coupling unification considerations it does not play a decisive 
role in establishing the correct neutrino mass scale. 

We first note that Fn has to be at least of the order of 0.2 x 10 -9 for 
the model to successfully reach a lower bound on the heaviest light neutrino 
mass that is approximately 0.05 x 10~ 9 GeV [32]. This is due to the well- 
known fact that b-r unification happens rather naturally in supersymmetric 
theories at the GUT scale. We find this to still be true even after we properly 
incorporate the intermediate scales in the running of fermion masses in our 
numerical studies when the SUSY scale is low. Namely, the b and r massess 
that are of the order of 1 GeV defer from each other by not more than 25 % for 
the allowed values of tan {3 and hence lead to a cancellation in the 33 element 
of the Md — M e matrix. (This sort of cancellation is required in order to get 
a large atmospheric neutrino mixing angle [US H3J H5] . ) Therefore Fn has to 
be about a factor of 4 bigger than one would naively expect in order to get 
a good fit of the overall neutrino scale. 

Clearly, a large enough Fjj prefers a low enough m barring some special 
cases when \fn(x)\ blows up. On the other hand, m is proportional to the 
GUT scale which is bounded from below by the experimental limits on proton 
decay. Therein lies the crux of the neutrino scale problem. Due to proton 
lifetime limits there exists a lower bound on m that implies an upper bound 
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on the attainable neutrino scale which, in general, tends to be too low to 
accommodate the experimental data on neutrino oscillations. 

To address this problem we use the current experimental bounds on the 
partial proton decay lifetimes to establish the phenomenologically viable 
lower limit on m. This allows us to establish an upper bound on Fu — and 
hence an upper bound on the neutrino scale — that the theory can generate 
for particular values of a, a, |A|, \r]\ and x through the type II contribution 
while being certain that proton is stable enough. We find that the following 
conservative limit on m applies 

m > 5 x 10 15 GeV|A|y^ ^4 

(6) 

For details of this derivation and notation we refer the reader to Appendix [A] 
As is taken to be a common value of the short distance enhancement factors 
of the relevant d = 6 proton decay operators between the GUT scale and Mz- 
These are calculable for a given mass spectrum of the theory. We take As 
to be 2.5 which overestimates the values obtained from running by 10-25%. 
Whenever we use Fn as defined in Eq. (|5]) with m explicitly replaced with 
the lower bound from Eq. ^ we refer to it as F^ ax . 

The last step allows us to systematically maximize Fu in order to find a 
region where it exceeds the desired value for a successful description of the 
neutrino data. Since _F™ ax ~ 1/ V% a small value of A is preferred. If and 
when such a region is found we can further check whether the gauge coupling 
unification takes place there. If it does, we generate the appropriate mass 
matrices by propagating the experimentally determined fermion masses from 
Mz to the GUT scale taking care of all the intermediate states, and use 
these values to see if a successful fit of fermion masses and mixing param- 
eters is indeed possible. This procedure guarantees a self-consistent check 
of the viability of the theory. (A similar line of approach has been taken in 
an SO(10) model with radiatively induced fermion mass hierarchy [IS] very 
recently [IT].) 

Before we optimize Fu, we first plot the lines of constant _p™ ax for one 
particular set of values of the relevant parameters (rj — X — a — a — 1, 
tan j3 = 50) in Fig. [T] to show how close to the right neutrino scale one 
usually gets. Note that we always run tan (3 (v) from the SUSY (M z ) scale 
to the GUT scale using the relevant equations to perform the numerical fit 
as well as to evaluate Fu- However, whenever we specify tan/3 throughout 
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this paper, we specify its value at the relevant SUSY scale unless stated 
otherwise. 

Lines of constant Fjf ax /10~ 9 shown in Fig. [I] correspond to 0.02 (solid 
line) and 0.002 (dashed line). Recall that these contours reflect the actual 
upper bounds on Fjj that are in accord with proton decay limits. Clearly, 
only in a very narrow region _p™ ax is within a factor of ten away from the 
preferred value. We can further ask if indeed unification takes place in that 
particular region and check whether it confirms our estimate for the mismatch 
between the required value for Fjj and the actual value for Fn. The dots 
in Fig. [T] mark places where successful one-loop gauge coupling unification 
takes place assuming that the GUT scale M GUT is given by the scale of the 
lightest proton decay mediating (X, Y) gauge bosons M(x,y)- We plot only 
the points that are in agreement with the experimental limits on proton decay 
lifetimes and yield a perturbative unification below the Planck scale. There 
are three different dot sizes; the largest dots represent unification scenario 
where F n /10~ 9 > 0.02, medium size dots correspond to 0.002 < F n /10~ 9 < 
0.02 and the smallest dots are for Fn/10~ 9 < 0.002. Here we note that not 
even one point corresponds to the case where Fn/10~ 9 reaches or exceeds 
0.2. In other words, we can be rather certain that no successful fermion mass 
fit with dominant type II scenario for neutrino masses can be found for these 
particular values of model parameters due to a too small neutrino mass scale. 

As far as the pure type I contributions are concerned we show in Fig. [T] 
that all the singular points of //, represented with squares, do not overlap 
with the viable unification points. In other words, whenever type I contribu- 
tion is potentially significant for light neutrino masses, unification does not 
happen. 



3.2 Constraints from unification and proton decay 

The exact procedure that we use to check for unification of gauge couplings 
for data in Fig. [T] is as follows. First, we fix the SUSY scale Msusy — 1 TeV. 
We further setr] = X = a = a = l and then vary x and gcur, where Qgut 
is the gauge coupling at the GUT scale. Once all these parameters are given 
we numerically determine the masses of the superheavy fields, including the 
gauge ones, in arbitrary units of m. Then we define the following coefficients 

B < = 2^ r/ ' ri = v~^T — m — ' o<?7<i, 7 

~ In MgutI Msusy 
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where b\ , i = 1,2, 3, are the usual one-loop coefficients of (3 functions of the 
7th threshold of the zth gauge coupling. Note that at this stage the only 
unknown is In Mgut/Msusy- Finally, we solve for it using the difference 
between the one-loop equations for the running of a± and a 2 from Msusy 
to M GUT and check that this value indeed generates unification and that the 
inferred gauge coupling at the GUT scale matches the input value for Qgut- 
This we do by requesting that (B 2 — B 3 )/(Bi — B 2 ) and oigut = QgutK.^) 
defer by less than 4% from the central values of (a 2 l — — a^ 1 ) 

and (c^ 1 — B 3 ln Mgut /Msusy / (Zn)) -1 , respectively. We introduce the 4% 
factors to reflect the fact that these are one-loop considerations for a single 
SUSY scale only and for central values of gauge couplings. We take a 3 = 
0.0895, a 2 = 0.0326 and a± = 0.0174 to be the input values at the SUSY scale 
of 1 TeV. These reflect the two-loop running effects from Mz to 1 TeV. We 
finally check that ccgut is perturbative, proton decay constraints are satisfied 
and GUT scale itself is below the Planck scale. To check the proton decay 
viability we look at p — > 7r°e + channel and the d = 6 operator contribution. 
To accurately evaluate the prediction for p — > 7i°e + we numerically find the 
relevant short distance coefficients, Asl and Asr, of the d = 6 operators 
(see Appendix [A| taking care of all intermediate scales using the results of 
Refs. [IHJ HHJ [50l [5TJ |52] and use Qguti M( X ,y) anc ^ ^(X',Y') as deduced from 
unification. 

In Figs. [2] and [3] we show the results of the fermion mass fit at the 
GUT scale when we actually run the charged fermion masses, the Cabibbo- 
Kobayashi-Maskawa (CKM) parameters and v from Mz to Mgut and tan /3 
from M S usy to M G ut for tan/3 = 30 and tan/3 = 50, respectively. When 
we run the charged fermion masses, the CKM parameters, v and tan/3 we 
use the knowledge of the full mass spectrum of the theory to properly in- 
clude all intermediate scales at the one-loop level. We do not run the light 
neutrino masses and the known angles of the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) mixing matrix. Again, the dots represent the points where 
gauge coupling unification takes place when r] = \ = a = a = l. The 
level of shading specifies the goodness of the fit, i.e., the range of values of 
X 2 , which exhibits clear dependence on tan/3. Since we already know that 
the correct mass scale for the neutrinos cannot be reached, we fit beside 
the charged fermion masses, CKM parameters, solar and atmospheric an- 
gles of the PMNS matrix, one particular ratio m 2 /m 3 of neutrino masses 
and then evaluate the mismatch between the mass of the heaviest light neu- 
trino m3 as inferred from the fit and lower bound on neutrino mass scale 
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that we take to be 0.05 eV. For the details of the fitting procedure we refer 
the reader to Appendix [Bj We find (0.05eV/m3) > 57 for numerical fits for 
which x 2 /d.o.f. < 10/13. This not only confirms our initial estimate that 
the theory fails to accommodate the relevant neutrino scale but gives us an 
accurate mismatch factor for the r] = \ = a = a = l case with low SUSY 
scale. 

3.3 The failure of the low energy supersymmetric case 

The main lessons to be taken from this example are the following. Firstly, 
Fu when combined with the limit on m as given in Eq. ^ provides direct 
means to accurately estimate whether a successful description of the light 
neutrino mass scale can be achieved. This should be the starting point of 
any particular scan of the parameter space of the theory. 

Secondly, m is just another parameter in the theory and as such must 
be allowed to vary as long as inequality in Eq. ^ is satisfied in order to 
cover all possible unification scenarios. It actually varies from 4 x 10 14 GeV 
to 4 x 10 17 GeV for unification points shown in Figs, [lj [2] and [3j Fixing m to 
a particular value would reproduce only a portion of the available parameter 
space. 

Thirdly, charged fermion masses vary drastically from point to point. The 
reason for this is that new (compared to the desert scenario) heavy states 
change the relevant renormalization group equations (RGEs). When we prop- 
agate the Yukawa couplings we neglect for simplicity the direct influence of 
GUT states threshold corrections, but take them consistently into account in 
the RGEs for the gauge couplings. This means that Yukawa couplings at the 
GUT scale differ from the ones obtained in the desert like MSSM scenario 
only due to a change in the values of the gauge couplings. We will comment 
on this later on. In Table [T] we present the ranges of the central values of 
masses of quarks and charged leptons after we extrapolate them to the GUT 
scale for the unification points presented in Fig. [T] and confront them with 
the central values when we assume MSSM-like scenario with tan/3 = 30, 
M GUT = 2 x 10 16 GeV and M SU sy = 1 TeV. This, in turn, affects the nu- 
merical fit which depends on tan/? for two reasons: it affects the RGEs as 
explained above and directly enters mass sum rules g-Q. 

Finally, it is clear that the unification consideration is crucial since unifica- 
tion happens in a rather narrow region of x. For example, all regions that are 
potentially viable in describing neutrino scale through the type I seesaw do 
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not overlap with the regions where couplings unify when 7] = X = a = a = l 
and are thus ruled out prior to any numerical fitting. 





m u 


m c 


m t 


m d 


m s 


m b 


m e 




m T 


min 


.00053 


.210 


86.7 


.0012 


.0213 


1.17 


.000341 


.0720 


1.28 


max 


.00060 


.237 


93.2 


.0014 


.0240 


1.29 


.000367 


.0775 


1.37 


MSSM 


.00050 


.198 


75.7 


.0011 


.0202 


1.07 


.000357 


.0754 


1.34 



Table 1: Range of central values of fermion masses in GeV units at the GUT 
scale when r] = X = a = a = l and tan (3 = 30 for different values of x 
where unification takes place. We also present the two-loop level running 
masses at the GUT scale in the MSSM-like setup when Msusy = 1 TeV and 
M GUT = 2 x 10 16 GeV. 

Clearly, the theory fails to describe the current data on fermion masses 
and their mixing angles with generic parameters. We thus resort to a different 
strategy that we already alluded to at the beginning of this section. Namely, 
we first maximize F™ x in terms of x, rj, A, a and a in order to make sure 
that _p™ ax is large enough for reasonable value of tan (3 to yield the correct 
neutrino mass scale. For that purpose we find suitable values of n, A, a 
and a by varying them in the following range: 0.03-7. Note that this fixes 
the value of m and thus sets the overall scale of the theory that is relevant 
for proton decay in terms of the scale that is relevant for the description of 
light neutrinos. Then we check whether unification of gauge couplings takes 
place. This we do by varying gcur within a given perturbative range to try to 
simultaneously satisfy all three RGEs of the gauge coupling constants while 
assuming that M S usy = 1 TeV. Again, once we know x, rj, A, a, a, m and 
9gut we have a full knowledge of the mass spectrum of the theory except 
Msusy- Finally, if and when unification works, we run the charged fermion 
masses, CKM parameters, tan/3 and v to the GUT scale and perform the 
numerical fit as described in detail in Appendix [B] to check the viability of 
the theory. 

Although the above analysis is an improvement of the existing studies 
[T31 dU [T5J [T6], we still do not find a satisfactory solution with a single, low 
SUSY scale. In fact, comes out short by a factor of 5 or more with respect 
to the minimal experimental value, primarily due to the fact that unification 
with a low SUSY scale does not allow for sufficiently small values of A. 
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4 Split SUSY case 



To reconcile the sufficiently small values of A with unification we split Msusy 
into two scales. The first scale stands for the common scale M1/2 of gauginos 
and higgsinos while the second one describes the common scale Mj of matter 
field superpartners. We now simply solve the three relevant RGEs for M1/2, 
Mj and gcuT to insure exact gauge coupling unification at the one-loop level. 
Within this scenario we find a multitude of solutions for various values of the 
parameters for which _p™ ax looks promising enough to attempt a numerical 
fitting of fermion masses and mixing parameters. As it turns out, these solu- 
tions are always of the split SUSY type with light gauginos and higgsinos in 
the 100 TeV range and superheavy squarks and sleptons that are close to the 
GUT scale (but still low enough to allow a supersymmetric treatment of the 
potential and of the mass spectrum). We accordingly perform the running of 
the Yukawa couplings from M z to M GUT taking care of the split supersym- 
metry effects as prescribed in Ref. [27] in the presence of intermediate scales 
for various values of tan f3 and execute a numerical fit of charged fermion 
masses, CKM parameters, two PMNS mixing angles and one particular ratio 
of light neutrino masses m 2 /m 3 . 

We find a solution that yields both successful unification and satisfactory 
X 2 which we explicitly spell out in Table [2] There we also specify our input 
values for charged fermions and CKM parameters at Mz- These Mz values 
are updates of results already presented in Ref. [53] to reflect the results of 
Ref. [53]. The numerical fit determines also leptonic angle sin#i3(= sf^ INS ), 
the only angle in the PMNS matrix that is yet to be determined. 

Our approach to the numerical fitting, as described in Appendix [Bj is 
statistical in nature. In view of that we cannot guarantee that our minimiza- 
tion procedure always finds a true global minimum of x 2 ■ 111 other words, 
there could still be some room for improvement as far as the numerical part 
of this study is concerned. Moreover, the complexity of the problem pre- 
vents us from scanning the whole available parameter space since we must 
always fix x, 77, A, a, a, m, gcur and tan (3 prior to numerical running and 
fitting of fermion masses and mixing parameters. Our main finding, though, 
remains: it is possible for the minimal renormalizable SO (10) to yield cor- 
rect neutrino mass scale while being in agreement with the gauge coupling 
unification paradigm as well as proton decay constraints. 

The minimal renormalizable SO (10) establishes a direct connection be- 
tween the neutrino scale and the GUT scale. In fact we have used this con- 
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1 11 

observable 


data at M z 


FIT 


RGE 


1 1 

pull 


m e (MeV) 




0.4866613 


0.42279 


0.42279 




nfy(GeV) 




0.10273 


0.08925 


0.08925 




m T (GeV) 




1.746 


1.534 


1.534 




m u (MeV) 




1.6 


0.56 


0.54 


+0.069 


m c (GeV) 




0.628 


0.218 


0.214 


+0.182 


mt(GeV) 




171.5 


72.0 


67.4 


+0.675 


m d {MeV) 




3.5 


0.24 


1.2 


-2.310 


m s (MeV) 




62 


27.2 


21.8 


+0.696 


mfj^Lrev ) 




2.89 


n 017 


u.yiu 


+U.U (\) 


Q CKM 
b 12 




0.2272 


0.2423 


0.2272 


+0.951 


a CKM 
23 




0.0422 


0.0478 


0.0474 


+0.199 


Q CKM 
5 13 




0.00399 


0.00447 


0.00448 


-0.015 


fiCKM 




0.995 


1.149 


0.995 


+0.598 


„PMNS 
b 12 






0.42 


0.55 


-1.091 


Q PMNS 
'-'23 






0.55 


0.69 


-0.764 


Q PMNS 
6 13 






0.103 






m 2 /m 3 






0.178 


0.180 


-0.104 



Table 2: Input (RGE) and output (FIT) parameters of the numerical fit at 
the GUT scale with x7d-cf. = 9.6/13, m 3 = 0.049 eV, m 2 = 0.0087 eV and 
mi = 0.0012eV. Unification takes place for \x\ = 0.109, arg(x)/ir = 0.52, 
a = 1.26788, a = 7, rj = 6.54112, A = 0.03, M 1/2 = 1.5 x 10 5 GeV, M~ f = 
9.0 x 10 13 GeV, Mqut = 5.8 x 10 15 GeV and g GUT = 1.3 (m = 6.5 x 10 13 
GeV). The input data at Mz are also given in the second column. 

nection as a guiding tool towards discovery of the viable parameter space. 
Since our solution requires a GUT scale that is rather low, any significant im- 
provement with respect to the current experimental limits on the p — ► 7r°e + 
lifetime would put serious constraint on this model. Of course, as already 
mentioned in the introduction, before reaching this conclusion, a two-loop 
check would be in order, due to the large sensitivity of the proton decay 
rate on the value of Mx,y and Mx\v- Be that as it may, we summarize 
in Table [3] the predictions of this model for the most relevant proton decay 
channels for the fit shown in Table [2j We derive them using the flavour 
structure dependence of the relevant d = 6 operators spelled out in Ref. |55j. 
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Partial 


mean life (10 33 years) 




p decay modes 


Fit: Table 


2 




Lifetime bounds: |56] 


Fraction (Ti/T) 


p — > 7r°e + 


1.6 


> 1.6 


44.3 % 


P — > 7T°H + 


70 


> .473 


1.0 % 


p -> K°e + 


442 


> .150 


0.2 % 


p -> irv + 


15 


> .120 


4.7 % 


p — > r/e + 


238 


> .313 


0.3 % 


p — > r//i + 


238 


> .313 


0.5 % 


p — > 7T°Z/ 


1.5 


> .025 


49.0 % 



Table 3: Predictions for partial lifetimes for the most significant p decay 
modes for the fit shown in Table [2] and current bounds [36] expressed in 
10 33 years units. For the matrix element we take a = 0.009 GeV 3 [57] . The 
predicted associated branching ratios are also shown. 



For completeness we explicitly provide the mixing matrices W e and V q , 
tan/3, and v at the GUT scale as well as Asl and Asr coefficients for the 
successful fit defined in Table [2] W e , V q , Asl and Asr are defined in Ap- 
pendix [Aj 

/ -0.9989 + 0.0260 i -0.0288 + 0.0227z +0.0136 - 0.0023 i\ 
W e = -0.0132 -0.0348 i +0.8338 + 0.5478 i -0.0559 - 0.0156 i , (8) 
\ -0.0096 -0.0080 % -0.0323 - 0.0486 i -0.7749 - 0.6293 i J 



I +0.8936 + 0.3778 i -0.2298 - 0.0771 i -0.0011 + 0.0043 i\ 
V q = -0.2386 - 0.0412 i -0.9648 - 0.0909 i -0.0438 + 0.0191 i\ (9) 
\ +0.0106 + 0.0002i +0.0440 + 0.0159 i -0.9828 + 0.1785 i J 

tan/? = 34.6, v = 147.9 GeV, A SL = 2.0, A SR = 2.2. (10) 

One comment is needed before we end this section. As we said before, the 
Yukawa couplings were run without taking into account the corrections due 
to new couplings appearing above the heavy thresholds. For small enough 
values of these new couplings, the effect is negligible due also to the small 
mass ratio between heavy thresholds and M GUT . But if the couplings are 
large enough (as they are here for example r\ and a), one may worry that 
they would completely change the Mqvt values of the light fermion masses 
and mixings that we then want to fit. We estimated this effect and found 
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out that it cannot dominate the running. There are more arguments in favor 
of our claim. First of all, the corrections to the 1-loop RGE due to new 
interactions above heavy thresholds is due to the wave-function renormal- 
ization of the light Higgs fields. This means that it is approximately flavor 
independent for short running and thus have little impact on the quantities 
we are fitting. Secondly, most of the lighter fields among the heavy ones are 
coming from 210 h- But the dangerous couplings can come from the opera- 
tors 210^126^126^ (large rj coupling) or 210^10^126^ (large a coupling). 
So the projections of fields coming from 126# or 126// are typically small in 
the light eigenvector directions, so the effective couplings are much smaller 
than the original SO(10) ones, and thus the effect becomes negligible. The 
same can be said for the 126 h contribution (a/77 suppressed) in the MSSM 
Higgs H U) d directions, which again causes a small coupling of these two fields 
to the weak triplet living in 210#. Finally, the contribution of the operators 
between SM singlets and the Higgs bilinears are small either because the 
components with lighter singlets have small couplings or components with 
larger couplings are heavy enough. 

5 Conclusions 

Of course we have not shown that the above solution is unique. In fact it 
is still possible that low energy supersymmetry is allowed by this minimal 
renormalizable SO(10): the unsuccessful fit to fermion masses and mixing 
parameters can be corrected by loops due to soft SUSY breaking terms. 
These same terms must be such to also reduce enough the typically large 
d = 5 proton decay operators. What we found in this paper is one very 
simple among such solutions: split supersymmetry with a well determined 
scale of 10 13 -10 14 GeV. Such a scale is large enough to sufficiently change the 
value of Qgut in order to allow for a small enough value of A which, in turn, 
yields the correct neutrino mass scale. 

Most predictions of this model are common to split supersymmetric mod- 
els: no sfermions at the TeV scale, the LSP a natural dark matter candidate, 
long gluino lifetime and the absence of flavour changing effects as well as 
d = 5 proton decay operators. What is further predicted in our model is 
the exact form and magnitude of the d = 6 proton decay operators: we 
know both the proton lifetime, together with all the branching ratios. Such 
a GUT is thus an example of a predictive theory of proton decay, similar to 
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non-supersymmetric theories, where the flavour structure alone (and not un- 
known supersymmetry breaking parameters) governs the proton decay rates 
|58j . On top of that we have a good description of all the fermion masses 
and mixings in the SM, connecting thus neutrinos with charged leptons. As 
a bonus of this fit we find out that the yet unknown parameter s^ INS turns 
out to be relatively large: 0.10. This cannot be treated as a prediction of the 
model, though. In fact, it is possible that other solutions exist. The same 
can be said for some of the other output numbers we have quoted above: for 
example, it is probably possible to find solutions with lower mass of some of 
the gauginos or higgsinos, if one allows them to be split. We confirm however 
that the low energy supersymmetric model with negligible corrections due to 
soft terms is not viable. 

We are aware of three weak points in the above analysis. 

First, some of the parameters are quite large (the biggest one is a = 7) so 
that perturbativity may be lost. Although (a/47r) 2 < 1, the large number of 
fields involved may make matters worse. Notice however that a large coupling 
is by itself not necessarily dangerous: to prove that the non-perturbative 
regime has been reached one would need to calculate the relevant processes 
(decay widths, cross sections) at the one loop level and compare them with 
the tree order expressions. Obviously such a huge calculation is beyond the 
scope of this paper. 

Second, as we already mentioned, we have not fully consistently run 
the Yukawa couplings; above the heavy thresholds only the gauge couplings 
RGEs were properly modified. Although we have estimated that the Yukawa 
couplings of our solution with split supersymmetry after running from Mz to 
Mqut do not get sizable corrections from the new operators above the heavy 
thresholds, it would be desirable to include these effects from the beginning. 
This would be needed among others to definitely prove that the low SUSY 
scenario is ruled out in this context. It is reassuring however, that in spite 
of the large couplings involved the actual corrections turn out to be small. 
This is a signal that the theory may still be in a perturbative regime, and 
that the first weak point mentioned above may not be crucial after all. 

Third, the split supersymmetry scale is pretty large and exceeds the 
usual cosmological bounds [59, 60J, although it is partially alleviated by the 
large gluino scale Mi/ 2 - Such a situation is similar to the case of radiative 
induced seesaw scale proposed in [61J: there also the split SUSY scale My 2 
was predicted to be large, close to the GUT scale, although the reason was 
different (to prevent right-handed neutrinos from becoming too light). Here 
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as there the solution to this problem could be a heavy enough gluino or a 
late inflation, so that gluinos are not produced later. This and the previous 
issues deserve a separate study, which we postpone for the future. 
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A Proton decay 

Here we derive a constraint on the mass spectrum of our model due to the 
current experimental bounds on the partial proton decay lifetimes. The idea 
is to replace the only relevant mass scale m in the model with the mass scale 
that can be inferred from experimental data. For the moment we assume 
that the d = 5 proton decay operators are subdominant to the d = 6 ones. 
If that is the case, the most stringent limit comes from the p — ► 7r°e + mode 
for which the current limit reads r(p —>■ 7r°e + ) > 1.6 x 10 33 years |54j . The 
model, on the other hand, yields the following width for that mode: 



^j- 2 A\ \a\\l + D + F) 2 [A* 8R ItfCWj) 11 + kl{V}f\V q Wir)\ 2 
+Al L kl\(W e ) n + (V q ) 11 (W e Vj) 11 \ 2 ] 1 (11) 



where W\ = U\U e , V q = U ] u U d and fc 1(2 ) = yJ2ira GUT /M(x,Y) {(x>,y>))- a sl(R) 
give a leading-log renormalization of the relevant operators from the GUT 
scale to M z . The QCD running below M z is captured by the coefficient A L . 
We note that the mass of proton decay mediating gauge boson M(x,y) in the 
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minimal SO (10) is given by the following expressions [2T] 



(XX) 



m- 



2x 2 + x - 1 


2 

+ 2 


x — 1 





2x(2x 2 + x - 1) 



(x 
rt 



(12) 



To progress we assume that A SL « = Ag, £7^ 
Mfx r ,Y>)- I n that case we infer the following conservative lower limit on m: 



t/e and M(x,y) 



m > 5 x 10 15 GeV|A 




2x 2 



x 



1 



X 



2x(2x 2 



x 



X 



I) 1 



-1/2 



To generate this result we use m p = 938.3 MeV, D = 0.81, F 



(13) 
0.44, U = 



139 MeV, A L = 1.25, \(V q ) u \ = 0.97377 and a = 0.009 GeV 3 [57]. 

Although the bound in Eq. (13) has been obtained after some assump- 
tions, we have checked at the end of the fit that the general expression for T 
in Eq. (11 ) is smaller than the experimental bound. 



B Fitting procedure in MSGUT 

Our approach to the numerical fitting of the relevant fermion masses and 
mixing parameters is straightforward. To perform the fit we first construct 
the minimization function of the input parameters, which we then numeri- 
cally minimize using a so-called simulated annealing procedure [62]. Since all 
the mass matrices depend on two complex symmetric Yukawa matrices, i.e., 
Y10 and Yj^, we have 15 real parameters at our disposal to fit all the fermion 
masses and mixing parameters. In practice, we trade Y w and for the 
down quark and charged lepton mass matrices to arrive at the following sum 
rule [33] : 

V q T M u V q = ^ tan/5 [(1 + £(x)) M d - ^{x)Wj M e w} . (14) 



N, 



d 



Here we work in the mass eigenstate basis of down quarks. We implement 
congruent transformation with the following convention: 

M x = UjM x U x , x = u,d, e, n, 

where M x are diagonal, real and positive. The input parameters we vary 
in the fit are the three down quark masses in Md and three angles and 
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six phases that parametrize W e = U\Ud- The charged leptons are basi- 
cally taken out of the fitting procedure due to their small experimental 
error (we fix them to the central value at the GUT scale). Therefore we 
are left with 12 parameters altogether, 6 of which are phases. Together 
with the model parameters that determine the spectrum (A, r/, a, a, x, m, 



Qgut) they completely determine the r.h.s. of Eq. (14) and determine the up 



quark masses and the CKM mixing matrix, which is multiplied by five addi- 



tional phases V q = W u U d = P x VckmPi, where P x = diag(e^ 1 , e^ 2 , e^ 3 ) and 
P 2 = diag(e i/3l ,e iA , 1). (Similar expression holds also for the PMNS mixing 
matrix Vg = U\U n = KiVpmns^-) To obtain the neutrino masses and V^, 
we simply plug and M e in the neutrino mass formula given in Eq. Q. 
We now specify the pull, Xh f° r eacn output value pi in the following way: 



Xi 



Pi ~Pi 
fiPi 



where Pi are masses and angles which have been run to the GUT scale, fi 
are the percentages that specify the errors and are listed in Table |4j The 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Pi 




m c 


m t 




m s 


m b 


C CKM 


a CKM 
b 23 


a CKM 
*13 


fiCKM 


fi (%) 


30 


10 


10 


35 


35 


10 


7 


4 


1 


20 



Table 4: Fitted parameters of the charged fermions and their relative errors. 

minimization function is then a sum of squared pulls and consists of the 
quark and the neutral lepton parts: 

10 / _ ~ \ 2 

2 2,2 2 l Pi Pi \ 

x = x, + x„ x, = E [-j^) • 

Again, we fix the charged lepton masses at their central values at the GUT 
scale. Therefore their contribution to the \ 2 is zero. Neutrinos, on the other 
hand, enter the fit in the following way. We notice that the fitting procedure 
results in a hierarchical spectrum of the light neutrinos, so we choose to fit 
the ratio of the two heaviest neutrinos m 2 /m 3 , which does not receive large 
RGE corrections. We further relax the maximal neutrino mixing angles, since 
the RGE corrections typically reduce them at the GUT scale in the case of 
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hierarchical spectrum in both type I and type II case [EH EH ES] to get: 

f m 2 /m 3 -0.18 V , ( s^ MNS -0.5 \ 2 
Xe ~ { O02 J + { 0.15x0.5 ) (15) 

/ PMNS_r ) a\ 2 / PMNS _ q i 7 \ 2 

/*23 0^\ /fig 0^\ 9 (4MiV5_ 017 A 

V o.ixo.6 J v 0.2x0.17 y v 13 ; 

As for the input values for all charged fermion masses at M z , they are 
taken from [53] and updated where needed to reflect the results of Ref. [54£] 
as shown in Table [2j In addition, we use the CKM mixing angles and the 
CP phase as inferred from the fit for the Wolfenstein parameters as given 
in [51|. For the gauge couplings at M z we take a 3 (M z ) = 0.1176 ± 0.0020, 
a 2 (M z ) = 0.033816 ± 0.000027 and a x {M z ) = 0.016949 ± 0.000005. These 
are then run at the one-loop level to the correct GUT scale for a given tan (5 
where we perform the fit. 
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Figure 1: Lines of constant F^ ax /10~ 9 for rj = X = a = a = l and tan/3 = 
50 and the corresponding viable gauge coupling unification points. Solid 
(dashed) line corresponds to F^ ax /10 -9 of 0.02 (0.002). Dots of varying sizes 
mark x values for which successful one-loop gauge coupling unification takes 
place for a fixed SUSY scale of 1 TeV. Their size indicates the magnitude of 
F n /10~ 9 . Large squares (triangles) mark the singular points of // (///). 
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^ 2 /d.o.f., tan /? = 30 
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|x| 

Figure 2: Result of the numerical fit of charged fermion masses, 1712/1713, 
CKM parameters, and solar and atmospheric angles for tan (3 = 30. Dots 
correspond to successful gauge coupling unification at the one-loop level for 
central values of low-energy observables when i] = \ = a = a = l and 
Msusy = 1 TeV. Their shading describes the goodness of the fit expressed 
through the magnitude of the \ 2 function. Notice that we do not perform 
the fit in the white region without dots, because of the absence of unification 
and/or too fast proton decay. 
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Figure 3: Goodness of the fit of fermion masses and mixing angles for 
tan f3 = 50, when r] = \ = a = a = l and M S ijsy = 1 TeV. 



26 



